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Abstract 



Fully chaotic Hainiltonian systems possess an infinite number of classical solutions 
which are periodic, e.g. a trajectory "p" returns to its initial conditions after some 
, fixed time Tp. Our aim is to investigate the spectrum {ti,T2,...} of periods of 

^ I the periodic orbits. An explicit formula for the density p{t) = ^ 5{t — Tp) is 

'nJ" ' derived in terms of the eigenvalues of the classical evolution operator. The density 

^ ■ is naturally decomposed into a smooth part plus an interferent sum over oscillatory 

O ■ terms. The frequencies of the oscillatory terms are given by the imaginary part of the 

complex eigenvalues (Ruelle-PoUicott resonances). For large periods, corrections to 
the well-known exponential growth of the smooth part of the density are obtained. 
An alternative formula for p(r) in terms of the zeros and poles of the Ruelle zeta 
^ ■ function is also discussed. The results are illustrated with the geodesic motion in 

>• . biUiards of constant negative curvature. Connections with the statistical properties 

k> ! of the corresponding quantum eigenvalues, random matrix theory and discrete maps 

^ I are also considered. In particular, a random matrix conjecture is proposed for the 

eigenvalues of the classical evolution operator of chaotic billiards. 
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1 Introduction 

In chaotic Hamiltonian systems, the unstable classical periodic orbits form a set 
of measure zero among all the possible trajectories. However, as has been empha- 
sized many times, the periodic orbits are of great interest. In particular, they are 
very important in the study of the structure of the phase space dynamics and the 
transport properties of the motion. Another relevant aspect of these orbits is their 
temporal behavior. At a given energy, the periods Tp of all the periodic orbits p 
form a discrete set, {Tp} = {ri, r2, . . . , Tj, . . .}. Many properties of this sequence 
are of interest in, e.g., understanding the quantum mechanical behavior of chaotic 
systems. For example, semiclassical theories establish a link between the statistical 
properties of the quantum eigenvalues and those of the classical periods {Tp}. One is 
led to answer questions like how grows the number of periodic orbits with increasing 
period, or what are the correlations (if any) between the periods of different orbits. 
In this respect, as regards the first question, one of the main results in the field 
is the exponential growth of their number with increasing period [1]. Concerning 
the second, there is no definite answer, but semiclassical arguments based on ran- 
dom matrix theory (RMT) suggest that there exists correlations between the orbits 
(orbits with similar period repel each other) [2]. 

Our purpose is to further explore the properties of the spectrum of the periods 
{Tp} of the periodic orbits of fully chaotic systems. We restrict the index "p" to the 
primitive orbits only, i.e. repetitions of a given orbit are excluded. Our main result 
is an explicit formula that relates the density of periods of primitive orbits. 



P(^) =$^^(^-^p) 



to the eigenvalues of the classical evolution operator, the so-called Ruelle-Pollicott 
resonances. These resonances, which generically are defined by a set of complex 
numbers (denoted {7}), characterize the decay of correlations in the time evolution 
of phase space densities [3] , and provide important information about the transport 
properties of the system. For an introduction see [4,5]. 

The set of periods {Tp} can thus be explicitly related to another, more fundamen- 
tal, set, the Ruelle-Pollicott resonances {7}. In this way, the properties of the r^'s 
can directly be related to the properties of the 7's. The results allow, in particular, 
to make a systematic analysis of the structure of the spectrum of periods, focusing 
from the larger towards the smaller scales. 

Semiclassical theories a la Gutzwiller or Balian-Bloch [6, 7] describe the density 
of quantum states in terms of the periodic orbits of the corresponding classical 
system. Here, in turn, we establish a connection between the density of periods of 
the periodic orbits and the Ruelle-Pollicott resonances, i.e. the eigenvalues of the 
classical evolution operator. In this way, our results allow to establish a link between 



the eigenvalues of the quantum and classical evolution operators that, hopefully, will 
be useful to more clearly elucidate their properties and correspondences. 

Not much is known at present about the distribution in the complex plane of the 
eigenvalues of the classical evolution operator of chaotic Hamiltonian systems (see, 
e.g., [4]). Paraphrasing Alfredo Ozorio de Almeida [8], it is fair to say that though 
of great theoretical interest, the formula to be derived below merely relates our 
ignorance of the periodic orbit spectrum to the even more mysterious maze of the 
eigenvalues of the classical evolution operator. However, there is an increasing effort 
to understand the properties and physical interpretation of the latter, and their study 
is a central theme in several of the most interesting recent developments in classical 
chaotic systems, and of their quantum counterparts. For instance, in quantum 
systems the Ruelle-Pollicott resonances lurk behind many interesting effects. They 
provide the most simple explanation of the long-range non universal correlations 
observed in the quantum spectra of bounded Hamiltonian systems [9,10], and show 
up in experiments that measure the transmission through open microwave chaotic 
cavities [11]. The spectrum {7} is also a central issue in recent field theoretic 
approaches whose aim is to demonstrate the validity of RMT in chaotic systems 
[12], and appears in some mathematical models of quantum chaos, the Riemann 
zeros [13]. More recently, several studies have clarified the correspondence between 
the quantum and classical propagators for discrete maps in the presence of noise [14]. 

The starting point of our study is an expression of the trace of the evolution 
operator as a sum over the periodic orbits (Q- An inversion of that formula, based 
on the Mobius inversion technique, is implemented. Assuming that the spectrum 
of the evolution operator consists of isolated resonances, the inversion leads to a 
general formula for the density p(r) in terms of the eigenvalues {7} (SJ. It is also 
shown, in the same section, that a natural decomposition of the density emerges, 
where resonances located on the real axis determine the average or smooth behavior 
of the density, while oscillatory interferent terms arise from the complex resonances. 
An alternative and mathematically more accurate formula for p(r), based on the 
Ruelle zeta function (instead of the determinant of the evolution operator), is first 
presented in ^ It serves as a reference for the calculations of §31 and complements 
the results obtained. Both approaches are compared in ^ Two illustrative ex- 
amples are worked out in ^ The first one is the geodesic motion in a billiard of 
constant negative curvature. The spectrum of the evolution operator can be explic- 
itly computed in this case. This allows to write down a formula for the density of 
periods of the periodic orbits, thus illustrating the general approach of 21 The re- 
sults reproduce, in our formalism, those of Ref. [15]. The second example is based on 
the Riemann zeta function. Though the dynamical basis for this model is hypothet- 
ical, it is included here mainly to stress the existing analogies and similarities with 
known results in analytic number theory. Finally, ^contains some general remarks 
and conjectures concerning the spectrum of the evolution operator, inspired by the 



results of [jHland quantum chaos theory. Special emphasis is put in the connection 
with the statistical fluctuations of quantum eigenvalues and random matrix theory. 
We also show that important qualitative differences exist between the spectrum of 
the classical evolution operator of smooth flows and that of discrete maps. 

2 The density of periods of periodic orbits: Ru- 
elle zeta 

In this section the aim is to derive an explicit formula for the density of periods of 
the primitive periodic orbits, Eq.(Q). The density will be expressed in terms of the 
zeros and poles, located in the complex plane, of a particular function, the Ruelle 
zeta function. To simplify, we will from now on make reference to the zeros and 
poles of the Ruelle zeta as its "singularities". The calculations presented in this 
section serve as a basis for those of |gl where the density p(r) will be expressed in 
terms of the eigenvalues of the evolution operator. 
The starting point is the function 

oo 

Pir) = Y.Y.^P^(^-'^p)- (2) 

p r=l 

The index r accounts for the repetitions (or multiple traversals) of a given primitive 
periodic orbit "p" . The function P{t) is naturally associated with a zeta function. 
To see this, we reproduce a well known derivation [4], which consists in including 
inside the double sum (j21) a factor exp[6(r — rTp)], whose value is one because of the 
delta function {b is a real positive constant). Then, expressing the delta function 
as 6{t — rTp) = (27r)~^ J_ exp[z^(r — rTp)]dC,, and making the change of variables 
^ = —iu, P{t) takes the form 

br °° rioo 

p J,— I J ioo 

For b positive and sufficiently large, the sum is convergent, and can therefore be 
interchanged with the integral. Doing that, and making the additional change of 
variables s = u + b, we now get. 

Finally, using the expansion log(l — x) = — Xl^i^"/"^' ^^ obtain the relation, 

nr) = -— / dse^^— log Zr{s), (3) 



where Zr{s) is the topological or Ruelle zeta function [16,17], 

ZR{s)=l[{l-e—^y' , (4) 

p 

for Re(s) large. 

Equation allows to make a connection between the sum (j21) and the analytic 
properties of Zji{s). In order to proceed we thus need some information about 
the latter. In fact, for certain classes of hyperbohc systems, it can be shown that 
Zr[s) is analytic for Re(s) > /itop, where /itop is the topological entropy, and has a 
meromorphic extension to the whole complex plane [17]. Restricting our analysis to 
this case, from Q an explicit and simple formula follows (with h > htop, and closing 
the path of integration towards the left part of the complex plane). It expresses 
P{t) in terms of the location in the complex plane, denoted rj, of the singularities 
of Zr{s), 

P{r) = J2g,e^^ . (5) 

V 

In this equation the integer g^j is the multiplicity, and is positive for poles and 
negative for zeros. The index rj has a double significance: it serves as an index to 
enumerate the singularities, and also denotes their location in the complex plane, at 

3 = 1]. 

Two alternative and distinct formulas for P(r) are thus available, one as a sum 
over the periodic orbits, Eq.©, the other as a sum over the singularities of Zji{s), 
Eq. (0) . It is convenient to rewrite Eq. Q in terms of the density, using the definition 

CO 

r=l 

From now on, the idea of the computation is the following. If we manage to invert 
equation (jHl), and express the density p{t) in terms of P{t), we are done, because 
Eq.(jni) can then be used for P{t), and the final output would be an expression of 
p{t) in terms of the singularities t] of the function Zii{s). 

Inversion problems have many important physical applications, but are in general 
difficult to solve. In our case, the inversion of Eq.® will be based on the Mobius 
inversion formula. This is a technique developed in the nineteenth century, that 
has been extensively exploited in number theory [18]. More recently, it has found 
concrete applications in physics, like for example for computing the phonon density 
of states from experimental measurements of the specific heat of solids [19]. 

The inversion proceeds as follows. Consider the sum 

S, = ± ^P(r/,„) . (7) 

^ — ^ m 

m=l 



where /i(m) is the Mobius function [18]. This is a number-theoretic function, whose 
properties are based on the prime decomposition of the integer m. It is defined as, 

{1 ii m = 1 

(—1)'^ if m is a product of k distinct primes 
ii m has one or more repeated prime factors , 

(and thus ^{m) = 1, —1, —1, 0, —1, +1, ... for m = 1, 2, 3, 4, 5,6,..., a quite erratic 
function). If, in Eg.©, we use for P the series defined by the r.h.s. of Eq.®, we 
obtain 
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•m,r=l ^ ' n=l m/n 

where the sum Ylm/n '^uns over the divisors m of n. The last equality is the key point 
of the inversion technique. It follows from a remarkable property of the Mobius 
function, namely J2m/nPi^ ~ ^n,i- Finally, combining the last equation with 
Eq.(I7j), we obtain an equation for the density. 
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The sum over m has been truncated at a value equal to the integer part of r JTmin-, 
where x^m is the period of the shortest periodic orbit of the system, 

"^c [''"/ Tmin\ ■ 

This truncation is a consequence of the fact that P(r) = for r < Tmin (cf Eq.©), 
and therefore Pir/m) = for m > uic in Eq.(jH)). 

Equation (jSI) is our first result. It expresses the density of periods of the primitive 
periodic orbits in terms of the singularities of the Ruelle zeta function. 

The right hand side of Eq.®, should, in principle, reproduce a series of delta 
peaks located, according to the definition of p{t), at the periods of the primitive 
periodic orbits. In order to better display the structure of Eq.®, and see the 
correspondence with the series of delta peaks, it is useful to consider a decomposition 
of the density in two parts, 

P{r) = Xr) + p{r) • (9) 

This decomposition is associated to a classification of the resonances into two distinct 
sets, each set contributing to one of the two terms in the r.h.s. of Eq.ijHj). Because 
P{t) is a real function, the Ruelle zeta satisfies a simple functional equation (which 
follows from Eq.®), 



(the star denotes complex conjugate). A singularity tj of Zr{s) is therefore either 
simple and real, or complex and comes in conjugate pairs symmetric with respect 
to the real axis. This property is at the origin of the decomposition Q, with the 
real singularities contributing to "pir), and the complex ones to ]}{t). If the location 
of the resonances in the complex plane is written as 

r] = qr,±itn , 

with q^ real and t^ real positive, then the two contributions to the density can be 
expressed as, 

PM = ;E^E9.^'""". (1°) 
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m=l rjGSR 



and 



P(^) = lll^ E 9, e"^-/- cos(t,r/m) , (11) 

m=l r], trj>0 

where in Eq. fjlO|) the sum is made over the real singularities of Zji{s), whereas in 
Eq. ()ll|) it is made over the complex ones located in the upper half part of the 
complex plane. 

The first contribution, p(r), is given by a sum of real exponential terms. It has 
therefore a smooth dependence on r, and describes the average properties of the 
density. It reproduces the behavior of the singular sum (^ when p(r) is smoothed 
on a scale which is large compared to the average spacing between delta peaks. 
To reproduce the average part of the density of periods we thus need to know the 
location of the real singularities. It is known, in particular, that in hyperbolic 
systems its rightmost real singularity, which controls the asymptotic growth when 
r — i> cx), is a simple pole at s = h^op [1], 

Vo = qo = htop ; 5-^0 = 1 . (12) 

The presence of this pole implies, keeping the term m = 1 in Eq. ()lU|l . 

Mr) ™ -^ . (13) 

We thus recover, to leading order, the well known exponential growth of periodic 
orbits in chaotic systems, with the typical rate given by the inverse of the topological 
entropy [1]. However, the inversion procedure employed here allows to go beyond 
that result, and obtain subdominant corrections to the average growth of the density. 
It is remarkable, indeed, that the contribution of the pole at s = htop gives also, 
from Eq. (jlUj) . a series of exponentially large corrections to the asymptotic behavior, 
whit signs depending on the Mobius function. 



1 ^ aim) f, ,^ e'^^^p^ ^h.^-.rl'i ^hto^r/^ 
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T ^-^ m r 2r 3r 

m=l 



The first corrections lower the density, while the first positive one occurs for m = 6. 

Besides the pole at /itop, other real singularities oi Zji{s), with real part g^ < htop, 
add further sub-leading corrections to the density of periods. Those located in the 
negative part of the real axis contribute with exponentially small corrections. We 
are not aware of any generic result about the location in the complex plane of the 
additional real singularities, and how their contributions compare to the corrections 
that come from 779. 

The remaining contribution to the density, p{t), behaves quite differently. This 
term is responsible for the discrete nature of the spectrum of periods of the periodic 
orbits. On top of the average behavior p, each complex singularity of Zr{s) adds 
an oscillatory term to the density of periods. It is the interference of the oscillatory 
contributions of all the complex singularities that, formally, reproduces the delta- 
peak structure of Eq.(|TI). 

The amplitude of each oscillatory term is given by the exponential of the real 
part g^ of the singularity. In contrast, the inverse of the imaginary part, 27rmt~^, is 
the period of the oscillation. Singularities with the smaller imaginary part describe 
long range fluctuations with respect to the smooth behavior of the density, on scales 
which may be much larger than the typical time that separates two neighboring 
orbits (terms with m > 1 are of longer range, but their weight is of lower order in 
the limit r -^ 00). As singularities with increasing imaginary part are included, 
details of p{t) on smaller scales are resolved. For instance, to distinguish individual 
peaks of p{t) located around r requires complex singularities with imaginary part 
of the order of the average density of periods, t^ ~ e'''°p'^/r. 

Eq.® provides therefore a harmonic decomposition of the density of periods, 
were the frequency of each sinusoidal wave is given by a fraction of the imaginary part 
of a complex singularity. Since arbitrary high frequencies are needed to reproduce a 
delta peak, Eq.® also implies that, generically, singularities with arbitrarily large 
imaginary part should exist. But we have no additional information about their 
distribution in the complex plane. 

Eq. (jHJ can be integrated with respect to r to obtain a harmonic formula for the 
cumulative distribution, or counting function of the periods. This function is defined 
as the number of primitive periodic orbits whose period is smaller than r, 

N{r) = [ p{x) dx = yQ{T-Tj,) , (15) 



where B is Heaviside's step function. The result of the integration is, 

yu(rri) 
m 

m=l rj 

where Ei is the exponential integral function. Similarly to Eq.(jH)), this expression 






is, formally, exact. The same decomposition and general remarks as for the density 
apply to this function. 

3 The trace of the evolution operator 

In the previous section, we described the periodic orbit density and the counting 
function in terms of the singularities of the Ruelle zeta function. In this and the 
next section we will introduce an alternative description of the density, based on 
the eigenvalues of the classical evolution operator. Our motivations for doing this 
are the following. On the one hand, an explicit relation between the periodic orbits 
and the eigenvalues of the evolution operator is of clear theoretical interest. It helps 
in understanding the properties of both sets, as well as their correspondences. On 
the other hand, the connection is a central issue in the study of the quantum and 
classical behavior of chaotic systems and allows, via semiclassical techniques, to 
relate the quantum and classical eigenvalues. 

Before presenting in ^the derivation of the alternative description of the density 
of periods, we need to previously introduce some definitions and basic equations 
related to the classical evolution operator. 

Consider a classical system whose dynamical state is defined by the phase-space 
coordinate x. After a time r, the point x evolves into a new state y = frix). The 
kernel of the evolution operator of the deterministic motion is defined by a Dirac 
distribution, 

CAy,x)=5iy-Mx)). (17) 

We shall assume throughout the paper that /r(x) describes the conservative smooth 
flow of a closed Hamiltonian system. Open systems are not treated here, but may 
be considered as well within the same formalism [4]. 
The trace of the evolution operator is defined as 



R(t) = / dx Cr{x,x) , (18) 

where x is integrated over the phase space. The function Cr{x, x) is the conditional 
probability density for the system to be at the point x by the time r if the initial 
state was at the same point. R{t) is therefore proportional to the classical return 
probability at time r. The trajectory starting at x and coming back to the same 
point after a time r defines a closed loop in phase space. Any phase-space closed 
loop defines a cyclic motion, because the system returns to its initial conditions. 
For a fully chaotic dynamics, an explicit expression of R{t) in terms of the periodic 
orbits of the system was obtained in Ref. [20], 

R(r) = YY-—-^ S{T-rTp). (19) 

^ ^ 4-erldet(M;-l)| ^ '^ ^ ^ 



The sum is made over the primitive periodic orbits p and over their repetitions r. 
Each orbit has a period Tp and a monodromy matrix Mp. The latter describes the 
stabihty of the orbit. The factor | det(Mp — 1)|~^ is related to the overlap between 
an initial cloud, centered initially around the orbit p, and its iterate, after a time 
r = TTp, in a linear approximation. The function R{t) has therefore a peak at the 
period of each primitive periodic orbit, or at one of its repetitions, with a weight 
inversely proportional to its stability. 

Notice that Eq. ()19|) is similar to Eq.Q, with the important difference of the 
stability factor in the denominator in the former. In spite of this difference, similar 
steps as those who led in the previous section from Eq.© to Eq.Q can be followed 
for R{t). They lead to a connection between the trace of the evolution operator and 
a zeta function. We will not repeat them here, but only mention that, for reasons 
to be understood below, the sign in the change of variables between ^ and uj is the 
opposite here. The result, analogous to Eq.Q, is, 

1 i-a+ioo r, 

Rir) = -— / ds e-^-- log Z{s) , (20) 

where a is a real constant < 0. The function Z{s) is called the spectral determinant, 
or Smale zeta function, and is defined as [4,16], 



Z{s) = exp 



E 



p.r 



r\ det{M; - 1) 



(21) 



The complex variable s has units of inverse of time. Eq. ()20|) relates the trace of the 
evolution operator to the analytic structure of the function Z{s). Under some general 
conditions valid for a certain class of hyperbolic systems (and analogous to those 
assumed in the previous section for Zji{s)), the Smale zeta function is generically an 
entire function [21,22] (i.e., analytic at all finite points in the complex plane). More 
complicated analytic structures of Z{s) may arise, like for example branch cuts that 
lead to power law decays in intermittent systems [9] . We will ignore here these other 
possibilities. Assuming, therefore, that Z{s) is entire, from Eq.()2U|) the trace of the 
evolution operator can be expressed as [20], 

Rir) = J29,e-^^ , (22) 

7 

where the 7's are the zeros of Z{s). The complex set of points {7} define the 
(complex and discrete) spectrum of the evolution operator. In this equation and in 
the rest of the paper, the index 7 has a double significance: it serves as an index to 
enumerate the zeros, and also denotes their location in the complex plane, at s = 7. 
The spectrum {7}, given by the zeros of the entire function Z{s), characterizes 
the relaxation towards equilibrium of classical statistical ensembles [3] . The complex 

10 



zeros are usually called Perron-Frobenius or Ruelle-Pollicott resonances. We shall, 
however, often refer to the whole set of zeros (real and complex) as resonances. 

The positive index g-y in ()22|) is the multiplicity of the resonance. The physically 
relevant modes being decaying ones, the corresponding resonances lie in the positive 
half plane Re(7) > (this justifies the choice of sign in the derivation of Eq.(j2nj). It 
is in general a difficult problem to determine analytically the spectrum 7 for a partic- 
ular system. However, in Hamiltonian systems, where the energy is conserved, there 
exists a well defined long-time equilibrium state, given by the invariant measure on 
the energy shell with the microcanonical weight. The existence of this equilibrium 
state is manifested in the analytic properties of Z{s) by the presence of a simple 
"ergodic" zero [5], located at 

70 = , 9-yo = l. (23) 

The ergodic zero, located at the origin, corresponds to the unique invariant mea- 
sure (in the sense of statistical ensembles) in fully chaotic systems. From Eq. (j22j) 
it follows that this very general property implies that liniT-^oo R{t) — > 1, which ex- 
presses the equiprobability to find the system at any phase space point. This fixes 
the normalization of the trace of the evolution operator, or "return probability" (the 
volume is set to one). Other resonances, with Re(7) > 0, are associated to decaying 
modes that describe the process of relaxation of an initial cloud towards equilibrium. 

4 The density of periods of periodic orbits: Ruelle- 
Pollicott resonances 

The two expressions of the trace introduced in the previous section, Eqs. (fTU|) and 
(j221); relate a sum over the eigenvalues of the evolution operator to a sum over all 
the periodic orbits of the system. Similarly to [JH we can exploit this connection to 
derive an explicit formula for the density of periods of primitive periodic orbits, p{t), 
but now expressed in terms of the eigenvalues of the evolution operator (instead of 
the singularities of the Ruelle zeta function). This formula will take again the form 
of a harmonic decomposition, i.e., the density will be expressed as an interferent 
sum over oscillatory terms. The frequency of oscillation of these terms will be given 
by the imaginary part of the Ruelle-Pollicott resonances. 

To compute p{t), now the starting point is Eq.(fT^. To express -R(r) in terms 
of the density of periods, the main obstacle are the stability factors, | det(MI' — 1)|. 
These are not only functions of the period of the orbits, but depend also on their 
stability. However, in fully hyperbolic Hamiltonian systems, where periodic orbits 
are dense in phase space, the stability or Lyapounov exponents of long orbits are 
the same for almost all orbits [23]. Therefore, for long orbits the stability factor 
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becomes only a function of the period, 

|det(M;-l)| ^-=^ |det(M;-l)|. (24) 

Within this approximation, and using the definition (^ of the density, the function 
R{t) may be rewritten as, 

oo ^ 

«W = -E,.|det(M;,^- 1)1 '>(-/'■)■ 

The factor | det(M^ ~ 1)1 is in fact a function of the variable rr (cf, for instance, 
Ref. [4]), and therefore | det(M^, — 1)| = | det(MT- — 1)| depends only on the period. 
The density then takes the form, 

oo ^ 

From now on, the situation is similar to that of ^ We should invert equation 
(J2S1), and express the density p(r) in terms of -R(r). Then, using Eq.(j221) for R{t), 
we will get an expression of p(r) in terms of the eigenvalues 7. It is convenient, to 
simplify the computation, to define a new function 

^, , |det(M^-l)| ^, , ^ I , , , 

G{r) = J ^-^ ^ R{t) = 5^ - p{T/r) . (26) 

Then, consider the sum, 

00 
5*2 = y^^fi{m)G{T/m)/r 






'm^ 



m=l 



where we again make use of the Mobius function fi{m). The sum 5*2 is evaluated in 
two different ways. First, using for G the series expansion in the r.h.s. of Eq. ()2(j|l . 
This gives, 

^ ^ y^ fi{m)p{T/rm) _ y^ ^ /i(m)p(r/n) _ 

m,r=l ^ ' n=\ m/n 

where the sum ^^ /„ is made over the divisors m of n. The last equality follows, as 
in ^ from the property J2m/n pi^ ~ ^'n,i- The second way to evaluate 5*2 is using 
for G the first equality in Eq. (j26|) . Combining both, we obtain, 

^ _ V- M!^ I det(M./^ - 1)1 R{T/m) . (27) 



T '■ — ' m 

m=l 
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As in ^ the sum over m has been truncated at a value equal to the integer part of 
r/Tmin) where r^m is the period of the shortest periodic orbit of the system. This 
truncation is a consequence of the fact that R{t) = for r < Tmin (cf Eq. flT^ ). and 
therefore R{T/m) = for m > rric in Eq. ()27|) . 

Eq.(j2Zj) gives an explicit connection between the density of periods and the trace 
of the evolution operator. To express the density in terms of the eigenvalues of the 
evolution operator, we simply use Eq.(j22I)- This gives, 

^ 1 y^ Mm) ^^^^^ _ ^)| y- ^-,r/m (28) 

m=l 7 

Equation (j28p is our main result. It provides an alternative formula for the 
density of periods of primitive periodic orbits. The structure of Eq. lj^Hj) is very 
similar to that of Eq.®, with two important differences: it includes the stability 
factors, and the sum is made now over the eigenvalues {7} of the classical evolution 
operator (given by the zeros of 2'(s)), instead of the singularities of ^_r(s). 

The right hand side of Eq.(j2Hl) should, in principle, also reproduce a series of delta 
peaks located, according to the definition of p(t), at the periods of the primitive 
periodic orbits. Since the Smale zeta function also satisfies the functional equation, 

(Z(.))* = Z(s*) , 

the density can be, as in ^ decomposed into two parts, p = J) + ]), where real 
eigenvalues of the evolution operator contribute to p, and complex symmetric ones 
to p. If the location of the resonances in the complex plane are written as 

with q^ and t^ real positive (or zero), then the smooth and oscillatory contributions 
to the density can be express as, 

_ ^ _ ^ Mmj ^^^(^^ _ ^)| ^ ^-,,.M , (29) 



r -^ m 

m=l 7g3?+ 



and 



^(^) = 7 E ^ I det(M./^ - 1)1 Y. 9i e-"^^/'" cos(t,r/m) , (30) 

m=l 7gC+ 

respectively (SfJ"^ denotes the positive real axis (including the origin), and C^ the 
upper right sector of the complex plane). 

How these results compare to Eqs. ljlOj) and (fTTjl ? Let's start the comparison with 
the average part. In the limit r — >■ cxd, the leading order contribution to Eq. ()29p 
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is given by the real zero with the smallest real part. That zero is known, it is 
the ergodic zero denoted 70 in the previous section, whose location at the origin is 
the only generic property known about the spectrum of eigenvalues of the evolution 
operator in hyperbolic systems. Keeping only the term q^ = and m = 1 in Eq. ()29|) . 

we obtain 

^.^00 |det(M.-l)| 
p{r) = . (31) 

Thus, in the limit r — > cxd, the leading order behavior of the average density of 
periods is proportional to the stability factor of the orbits. When compared to 
the leading order obtained in ^ from the leading pole of the Ruelle zeta function, 
Eq. (fT^ . this result implies, 

|det(M, -l)r=°"e'''°p" . (32) 

This correspondence is equivalent to the Hannay-Ozorio de Almeida sum rule, de- 
rived in [23,24] from a uniformity principle. It expresses the counterbalance between 
the exponential proliferation of the periodic orbits and the growth of their (positive) 
instability. It is also a consequence of Pesin's equality, that relates the topological 
entropy to the sum of the positive Lyapounov exponents [25] . Here it has an analytic 
significance, it expresses the correspondence between the pole at htop of the Ruelle 
zeta function and the ergodic zero at the origin of the Smale zeta function. 

The higher order terms m > 1 obtain from the ergodic zero in Eq.()29|) give for 
the average density a contribution equivalent to Eq. (fT^ . Besides the ergodic zero, 
other real resonances (g^ > 0,t-y = 0) of Z{s) add further sub-leading corrections 
to the density of periods. However, no generic result concerning their location is 
known. 

The remaining contribution to the density, p(t), gives the oscillatory terms. In 
this approach, based on the eigenvalues of the evolution operator, the amphtude 
of each oscillatory term is determined by a competition between the exponential 
growth of the stability factor | det(Mr/m — 1)| and the exponential decay of e""^^"^/™. 
Using the asymptotic approximation ()32|) for the determinant, general arguments 
indicate that the complex resonances with the smaller real part should be located 
on the strip < g-y < /itop, and therefore that the leading oscillatory contributions 
have individual amplitudes that grow exponentially in time [26] . The inverse of the 
complex part, 2'Kmt~^, is the period of the oscillation. Resonances with the smaller 
imaginary part describe long range fluctuations on top of the smooth behavior of the 
density, on scales that can be, depending on the value of t-^, much larger than the 
typical time that separates two neighboring orbits (terms with m> 1 are of longer 
range, but their weight is of lower order in the limit r -^ 00). As resonances with 
increasing imaginary part are included, details of p{t) on smaller scales are resolved. 
For instance, to distinguish individual peaks of the density located around a period r 



14 



requires complex resonances with imaginary part of the order of the average density 
of periods, t-^ 



JltopT 



Ir. 

Eq. ()28p therefore provides an ahernative harmonic decomposition of the density, 
were the frequency of each sinusoidal wave depends on the imaginary part of the 
position in the complex plane of the eigenvalues of the evolution operator (Ruelle- 
PoUicott resonances). Since arbitrary high frequencies are needed to reproduce 
a delta peak, Eq. (j28|) implies that, generically, resonances with arbitrarily large 
imaginary part should exist (similarly to the Ruelle zeta function). 

Eq.()28|) can be integrated with respect to r to obtain an alternative formula 
for the cumulative distribution, or counting function of the periods, Eq. (fT3|) . To 
compute the integral, an explicit form of the determinant is needed, that depends 
on the dimensionality of the system. However, asymptotically the approximate 
expression (jH^ can be used. This yields. 



N{ 



T 



rric / \ 



m=l 



'-top 



7): 



ml 



(33) 



which has to be compared to Eq.()16|). The same decomposition and general remarks 
as for the density apply for this function. The main contributions to the smooth 
part are obtained from the ergodic zero, 7 = 0, and oscillatory terms are given by 
the complex resonances. 

We have therefore obtained two alternative and distinct descriptions of the den- 
sity of periods of the periodic orbits. From a mathematical point of view, Eq. (jHJ is 
preferable to Eq. ()28|) . because in the derivation of the latter we have ignored some 
possible fluctuations of the stability factors that may occur at short times. However, 
from a physical point of view the latter description is clearly more interesting, since 
the Ruelle-Pollicott resonances, or eigenvalues of the classical evolution operator, 
have a clear and intrinsic physical content, directly related to the classical dynamics. 
In practice, it is interesting to exploit both approaches. This leads to consider in 
more detail the relationships between them. 

As we have shown, the ergodic zero of Z{s) and the pole of Zr{s) at /itop carry 
similar (though not exactly equivalent) information. The pole at /itop leads to the 
smooth contribution (|14|). whereas the ergodic zero generates a series which is iden- 
tical but with e'^*°p^/'" replaced by | det(MT-/m — 1)|- It is only in the asymptotic 
approximation ()32|) of the stability factor that their contributions coincide. It is in- 
teresting to explore the correspondence between the analytic structure of both func- 
tions in the latter approximation. The simplest procedure is to make, in Eq. (j2ip . 



the replacement 
takes the form, 

Z[s) ~ exp 



det(M; 



ito-pr 



This approximate Smale zeta function 



E 



p,r 



3(s-fi.top)rTp 



n[i 



_ e(s-h.op)r,] =z-^\h,,^-s) 



(34) 
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where we have used the expansion log(l — x) = — Yl^=i ^"^ 1^^ to compute the sum 
over the repetitions. In this approximation, the poles/zeros r] of Z^{s) are mapped 
into zeros/poles 7 of Z{s), located at 7 = /itop — ^- This correspondence is expected 
to hold for long times, where the approximation of the stability factor made before 
is valid. Using the location of the singularities of Zji{^s\ the mapping (jH^ should 
therefore accurately describe the location in the complex plane of the zeros of Zi^s) 
that are responsible for the asymptotic behavior of the density. Those are the zeros 
whose real part is small. It maps, in particular, the pole of Zjii^s) at /itop into the 
ergodic zero of Zi^s) at the origin. In contrast, resonances with real part far away 
from the origin need not be in correspondence. Moreover, the approximation ()34j) 
produces a function Zi^s) which is meromorphic, rather than entire, as assumed 
before. Globally, the mapping (j34j) is clearly false. 

5 Illustrative examples 

The aim now is to find systems were the analytic structure of Zi^s) and Z^i^s) can 
be computed, thus allowing to write-down explicitly the formula for the density 
of periods of the periodic orbits, and to compare and illustrate both approaches. 
Another interest of such a study is to gain some insight into the distribution of the 
resonances and singularities in the complex plane in concrete examples. The explicit 
computation of the analytic structure of the zeta functions is a quite difficult task 
in general. However, it is doable in some cases. Two examples are treated in detail. 
The first one is the geodesic motion on a two-dimensional manifold of constant 
negative curvature. The second, a mathematical model, is based on Riemann's zeta 
function. Both have their own peculiarities: the somewhat unphysical character 
of the motion in the first, the purely speculative dynamical interpretation in the 
second. In spite of them, they both provide a concrete illustration of the general 
results obtained previously. The geodesic motion on a space of negative curvature 
has long served as a paradigm of classical and quantum chaotic motion [6,27]. 

5.1 Geodesic flow on surfaces of constant negative curvature 

We are interested in the geodesic motion on a two-dimensional hyperbolic geometry. 
In particular, we consider billiards on surfaces of constant negative curvature, the 
so-called Hadamard-Gutzwiller model. We will not enter into a detailed description 
of the classical and quantum motion in such surfaces, and refer the reader to the 
excellent introductory articles [27]. For these models Gutzwiller's trace formula is 
exact [28]. This allows to express the periodic orbit length spectrum in terms of 
the quantum eigenvalues [15]. We will here re-derived the formula for the length 
spectrum as an illustration of the general formalism based on the eigenvalues of the 
evolution operator and singularities of the Ruelle zeta. 
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Consider a compact and closed surface of area A on the two-dimensional Poincare 
disk, constructed from a suitable bounded domain on which appropriate (periodic) 
boundary conditions have been defined. The corresponding classical geodesic mo- 
tion includes a set p of primitive periodic orbits (closed geodesies on the compact 
surface). To start with, consider the corresponding Selberg zeta function 



Zs{s) = l[l[{l-e-^'+''^'^) . (35) 

p n=0 

The first product is defined over the primitive periodic orbits p of length ip (in 
appropriate units). This function differs in its definition from the Smale and Ruelle 
zeta functions. Zs{s) is an entire function, with zeros located at [29] 

a) s = l; g = l 

h) s = 1/2 ± ipa] g = ga 

c) s = 0; g = A/{2n) + 1 

d) s = -k, A; = 1,2,...; ^ = (A; + l)^/(27r) , 



where Pa = \/Ea — 1/4, a = 1,2,.... The Ea are the (quantum) eigenvalues of 
the Laplace-Beltrami operator on the surface, and the Pa are the corresponding 
wavenumbers. The integer g denotes the multiplicity of the zeros. In case b) they 
depend on the degeneracy (^q > 1 of the eigenvalue Ea- Because of topological 
constraints. A/ {An) is a positive integer. 

The topological or Ruelle zeta function (JH) is easily obtained from Zs{s). It is 
given by [29], 

Z„W=n(l-e-"')-' = ^fj±j^. (36) 

The analytic structure of this meromorphic function directly follows from that of 
Zs{s) (see the left part of Fig.l), 

a) pole at s = 1; g-q = ^ 

b) poles at s = 1/2 ± ipa] gr, = ga 

c) pole at s = 0; (?,, = A/{2'k) 

d) zeros at s = —1/2 ± ipa, g^q = —ga 

e) pole at s = —1; (7^ = A/{2n) — 1 

f) poles at s = —k, k = 2,3, . . .; g^^ = A/{2-n). 

17 



The rightmost pole of this function is real and located at s = 1. This implies 
that htop = 1 in these systems. This pole is responsible for the leading asymptotic 
average growth of the number of orbits, and provides sub-leading corrections as 
well. There is an infinite number of other real poles with g^ < htop, which also 
provide sub-leading corrections to the average part. The complex zeros and poles, 
aligned here on two different vertical lines in the complex plane, contribute to the 
oscillatory part. There exist poles and zeros with arbitrarily large imaginary part, 
as generically required (cf Q. 

Taking into account separately the contributions of the real and complex singu- 
larities of Zii{s), from Eq.® a formula is obtained for the length spectrum of the 
periodic orbits on a compact and closed hyperbolic surface (with the speed set to 
one, r stands here for the length i of the orbits). 



with 
and 



p{r) = p{r) + p{t) , r = length of orbits, (37) 

(38) 



/^("^) ^r/m 



e 



r ■' — ' m 

m=l 



A p-T/m 



27r d - e-^/™) 



A ( \ 

pij) = -y^. sinh(r/2m) Y^ 5f„ cos(p„r/m) . (39) 



T -^ — ' m 



Concerning the complex singularities, as mentioned before Zr^s) has two "crit- 
ical" lines, one made of complex zeros (located at Re(s) = —1/2), the other of 
complex poles (located at Re(s) = 1/2). Their superposition produces the oscilla- 
tory part ]){t). It is remarkable that the frequencies of the harmonic decomposition 
of the periodic orbit density are directly related to the quantum wavenumbers, a 
consequence of the fact that Selberg's trace formula is exact. We will come back 
to this point later on. The integrated version of Eq. ()37p for the counting function, 
N{t), was computed and analyzed in [15]. For a specific hyperbolic billiard, we 
have numerically checked that the deviations with respect to the leading order be- 
havior N{t) ^ Ei(r) observed in the data were well explained by the sub-leading 
corrections m > 1 of the pole r] = htop = 1 in Eq. (fT^ . 

Notice that each of the terms in the series (jHHj) diverges in the limit t —>■ 0. This 
divergence is due to the splitting of the density into two parts, p and p, whereas the 
sum of the two terms is well behaved. 

Consider now the spectral determinant or Smale zeta function ()2H) . whose ze- 
ros are the eigenvalues of the evolution operator. Another peculiar feature of the 
Hadamard-Gutzwiller model is that all the periodic orbits have the same Lyapounov 
exponent (equal to one). Mp can be written as a two by two diagonal matrix whose 
diagonal elements are e^^p. Therefore, 

I det(M; - 1)1 = 2 [cosh{rip) - 1] . (40) 
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Expanding the inverse of the determinant as | det(Mp — 1)| ^ = YlT=i ^ ^ ^ ' ^^'^ 
computing the sum over the repetitions r, from (J2H) and (j4(jp we obtain 

oo 

Zis) = llIl[^---''-''''T ■ (41) 

p k=l 

Using the definition of the Ruelle and Selberg zeta functions, Eq. ^Tj) can be re- 
expressed as 

oo oo ykfi \ 

Further manipulations of the latter equation lead finally to 

oo 

Z{s) = l[Zs{k-s) . (43) 

fc=i 

Since Zs{s) is entire, it follows from the last expression that Z{s) is also an entire 
function. The analytic structure of Z{s) follows from that of Zs{s), with zeros at 
(cf the right part of Fig.l), 

a) s = 0; g-y = I 

b) s = k + l/2± ipa, k = 0,l,2,...;gy= g^ 

c) s = k, fc = l,2, ...; g^ = k{k + l)A/{A7T) + 2. 

This distribution of zeros satisfies the "generic" requirements concerning the spec- 
trum of the classical evolution operator of an hyperbolic system: (i) all zeros have 
1^6(7) > 0, (ii) there is a simple pole at the origin (ergodic zero), and (iii) there are 
complex symmetric zeros with arbitrarily large imaginary part. Using Eq. ()4()|l . the 
density ()28|) is now expressed in terms of the classical resonances as, 

p^r) = 1 f] ^ e-/- (1 - e--/'")' Y. 9, ^''^'"^ ■ (44) 

m=\ 7 

Using the locations and corresponding multiplicities of the zeros of Zi^s) given above, 
and separating the contributions of the real and complex zeros, it is easy to check 
that Eq.(jlH) is strictly equivalent to Eqs. fjHHj) and (jSH). For the geodesic flow on 
a constant negative curvature, the exact density is thus recovered from Eq.ipHj). 
without any error. The reason for this is that in the present model the Lyapounov 
exponents are constant, independent of the orbit. It follows that Eq. ()24|) is exact 
for any period, not only asymptotically (cf Eq. pUj) ). 
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Figure 1: Analytic structure close to the origin of the Ruelle zeta function (left) and 
of the spectral determinant (Ruelle-Pollicott resonances, right) for a compact and closed 
billiard on a surface of constant negative curvature. Crosses are poles, circles are zeros. 
Multiplicities and exact positions are given in the text. 



It is however interesting to remark that the two expressions for the density are 
obtained from functions with very different analytic structure. In particular, Zii{s) 
is meromorphic and has two "critical" lines, one made of an infinite number of 
zeros located at Re(s) = —1/2, the other of an infinite number of poles located 
at Re(s) = 1/2. In contrast, the spectral determinant is entire and, as complex 
zeros, has an infinite number of parallel replicas of the quantum spectrum located 
at Re(s) = k + 1/2, A; = 0, 1, 2, . . .. The correct result (jHIJjl for the oscillatory part 
of the density is only recovered when the whole set of complex resonances of Z{s) 
is taken into account. In contrast, the leading order asymptotic behavior r -^ oo 
of the density is controlled by the zero at the origin and the complex resonances 
located at Re(s) = 1/2. 

The Hadamard-Gutzwiller model offers also the opportunity to analyze the accu- 
racy of the approximation (jMj) of the Smale zeta Z{s), which holds asymptotically. 
It is obtained by keeping only the first term fc = 1 in Eq. ()42|) . In that approximation, 
the analytic structure of Z{s) consists now of a simple zero at the origin, complex 
zeros on the line Re(s) = 1/2, complex poles on the line Re(s) = 3/2, and zeros at 
s = k,k = 1,2,.... The analytic structure is thus well reproduced for Re(s) < 1 
(the zero at the origin and the lowest line of complex zeros located at Re(s) = 1/2). 
The remaining structure is wrong (the degeneracies of the remaining real zeros is 
wrong, at Re(s) = 3/2 it has a line of complex poles instead of zeros, the other lines 
of complex zeros are missing). 
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5.2 The Riemann zeta 

Our second example is taken from analytic number theory, which has inspired several 
developments in the theory of dynamical systems [30] . The results we are going to 
present are based on a dynamical interpretation of the Riemann zeros and of the 
prime numbers. This interpretation is by no means necessary, but it is a useful one 
because it introduces the appropriate physical framework into the discussion, and 
therefore facilitates the comparison with dynamical systems. We therefore revisit 
here some well known formulas in number theory, viewed from the perspective of the 
present theory that connects the periodic orbits to the eigenvalues of the evolution 
operator. 

The spectral interpretation of the Riemann zeros is based on the following iden- 
tification. The imaginary part of each of the Riemann zeros is thought to be an 
eigenvalue of a quantum system with a classically chaotic limit with no time-reversal 
symmetry. An analysis, based on a semiclassical interpretation of a formula for the 
density of the critical zeros [31], shows that the set of prime numbers has to be 
identified to the set of periodic orbits of the (unknown) classical dynamics. The 
analysis leads to the following mapping, 

primitive periodic orbits -^ prime numbers p 
period of the orbits Tp — > logp (45) 

stability I det(M; - 1)| ^ f = e"^°sP _ 

The last correspondence implies that /itop = 1 in this hypothetical dynamical system, 
that for simplicity we refer to as the Riemann dynamics. 

We don't know the classical Hamiltonian behind the Riemann dynamics, but 
the information contained in ()45|) concerning the periodic orbits is enough to write 
down the trace of the corresponding classical evolution operator, Eq.()19|). Using in 
the latter the correspondences (j^ . the trace is expressed as, 

oo 1 2 
p r=l 

where the sum runs over the prime numbers p. To write down R{t) in terms of the 
classical Ruelle-PoUicott resonances, we need to compute the corresponding spectral 
determinant, or Smale zeta function. The result, obtained from (|2H) and (|45p is [13], 

Z{s)=C\l-s), (47) 

where ({s) = Yl'^=i^~^ i^ ^^^ Riemann zeta function. Contrary to general ex- 
pectations in bounded hyperbolic systems [21], Z{s) is not entire for the Riemann 
dynamics, but meromorphic. The analytic structure of Z{s) follows from that of 
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a) zero at s = 0, g-y = I 

b) poles at s = 1/2 ±i ta, g^y = —ga 

c) poles at s = 1 + 2k, k = 1,2, . . ., g^ = —1. 

Assuming the Riemann hypothesis, the ta are real and define the position of the a'th 
zero of ({s) along the critical line, of multiplicity (7^ > 1. The pole of the Riemann 
zeta transforms into the ergodic zero of Z{s). The position of the complex Ruelle- 
PoUicott resonances coincide here with the complex zeros of the Riemann zeta, but 
these are not zeros of Z{s), but poles. Other real poles of Z{s) are generated by 
the so-called trivial zeros of the Riemann zeta, now located on the real positive 
axis. The global properties of the distribution of the singularities of Z{s) for the 
Riemann dynamics satisfy the general requirements of fully chaotic systems, with 
the important oddness related to the occurrence of some poles in place of zeros. 

From (j^ and the analytic structure of Z{s) (assuming ga = 1), the function 
R{t) may be written, 

Rir) = 1 - . ^ _2, - 2 e-^/2 J2 cositaT) . (48) 

The "anomalous" minus signs that appear in the r.h.s. of this equation refiect, 
again, the presence of poles in Z{s). The physical interpretation of these signs, and 
of the closely related minus sign that appears in front of the oscillatory part of the 
density of the Riemann zeros [31], is unclear for the moment, though an appealing 
possibility was suggested in [32]. 

The function R{t), here interpreted as the trace of the classical evolution oper- 
ator, is in fact a well known function in number theory, 

R{t = \ogx) = — — , 
ax 

where ?/'(x) = Xln<a;^(^) (^l'^) ^^ "^^ Von-Mangoldt function) [18]. By inverting 
Eq. ()46|) , as in ^ and using Eql48|) , a formula for the density of the logarithm of the 
prime numbers is obtained, given by Eq. ()28p with the appropriate identifications 
and replacements dictated by the table P^. The result coincides, up to a change 
of variable, with Riemann's formula. It can be integrated to obtain the counting 
function. These results stress the strong similarities that exist between number 
theory and the theory of dynamical systems. 

The Ruelle zeta function (0)) coincides with the Riemann zeta [30], Z^i^s) = 
Yip (1 — P~^)~ = C{s)- The approximation ()34|) is therefore exact in this case, since 
Eq.dnai is (cf Eq.dllD). 
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6 Concluding remarks 

Two different explicit formulas for the density p(r) of periods of the primitive pe- 
riodic orbits of fully chaotic classical systems have been obtained. Both provide a 
harmonic decomposition of p{t), where the complex zeros and poles (if any) of the 
corresponding zeta function are related to the elementary frequencies of the oscilla- 
tory terms, while the real ones contribute to the smooth part. In one case, the zeta 
function is the spectral determinant Z{s). It is assumed to be an entire function; 
its zeros, denoted 7, are the eigenvalues of the classical evolution operator (usually 
called Ruelle-Pollicott resonances). In the second formulation, the relevant function 
is the Ruelle zeta Zji{s). In contrast to the spectral determinant, this function has a 
meromorphic extension in the complex plane. The relation between both approaches 
was discussed in some detail. 

The zero of Z{s) located at the origin (or, alternatively, the pole of Zu{s) located 
at s = /itop) provides the leading average growth of the density of periods of periodic 
orbits, IJiT) ~ e^'°p^/r. We found exponentially large sub-dominant corrections 
to the leading term that arise from the same zero (or pole), and that are responsible 
for the main deviations observed numerically in billiards in a surface of constant 
negative curvature. 

The zero of Z{s) at the origin reflects a generic property of hyperbolic systems, 
the existence of an equilibrium distribution described by the micro canonical mea- 
sure. What about the rest of the spectrum? No generic statements are known, 
aside the fact that Re(7) > 0, that resonances with arbitrarily large imaginary part 
should exist, and that complex resonances come in symmetric pairs (cf Q. There 
are, however, some hints on what could probably be the generic structure, if any, of 
the low lying spectrum of Z{s) (e.g., resonances with the smaller real part), that we 
would like to discuss now. This low-lying part of the spectrum describes the long 
time dynamics of the system. For simplicity, from now on we restrict the discussion 
to the case of billiards. We have in mind "generic" systems with a discrete spectrum 
of the classical evolution operator (exponential decay). We also restrict to ballistic 
systems (i.e., bilhards whose shape produce a chaotic motion), and are excluding 
disordered systems. 

In ^ we saw that for billiards on constant negative curvature the low-lying 

spectrum consists of the ergodic zero at 70 = 0, and of an infinite number of zeros 

located on the line Re(s) = /itop/2. This structure also follows for other systems from 

semiclassical arguments. For a fully hyperbolic billiard, Gutzwiller's trace formula 

for the density of quantum eigenvalues takes the form, 

1 J- 

Pq(^) =y2^^^~ ^") ~ ^q(^) + ~ 5Z / "^ = cos(rtrp) . 

a ^7^^|det(M;-l)| 

As in ^5.11 Tp denotes here the length of the periodic orbits p and ta are the (eigen) 
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wavenuinbers; 'Pq{p) is the Weyl term, and we have neglected Maslov indices. Fourier 
inverting this formula with respect to t, and using the approximation (j32|) . we obtain 
the following formula for the density of periods, 

p(r) ^ + - e^*°p^/2 y- cos(t„r) . (49) 

r r ^-^ 



This is precisely the leading order density that is obtained from (j28|) assuming that 
Z{s) has a simple zero at the origin and complex zeros concentrated on the critical 
line Re(s) = /itop/2. Thus, the Hadamard-Gutzwiller model of [jSl as well as the 
generalization Eg . (j49j) suggest that the generic low-lying spectrum of the classical 
evolution operator of fully chaotic billiards consists of a simple zero at the origin 
plus an infinite sequence (extending to arbitrary large imaginary parts) of complex 
symmetric zeros located on the line Re(s) = /itop/2. A corresponding structure 
follows for Z^{s) by transforming zeros into poles located at /itop ~ 7- This first 
hypothesis about the location of the complex zeros of Z{s) is reminiscent to that of 
Riemann in number theory. 

For the negative curvature model and in the inverse formula (jm|) . the imaginary 
part of each Ruelle-PoUicott resonance located on the line Re(s) = /itop/2 coincides 
with a quantum wavenumber. This happens because of important non genericities 
of those systems: the corresponding semiclassical trace formula (the Selberg trace 
formula) is exact, and Maslov indices are zero. If in Eq. (j49p the Maslov indices are 
not neglected (and correction terms are taken into account), the connection with 
the quantum wavenumbers is generically lost. The simplest effect of these phases 
would be to produce a shuffling of the resonances, without moving them out of 
the line Re(s) = /itop/2, and without changing their statistical properties. Some 
arguments in favor of this will be given below. Concerning the distribution of the 
eigenvalues of the classical evolution operator, our second guess is therefore that 
asymptotically (i.e., for resonances located far from the real axis) the statistical 
properties of the imaginary part of the Ruelle-PoUicott resonances located on the 
critical line Re(s) = /itop/2 coincide with those of the corresponding quantum (eigen) 
wavenumbers, which are random matrix like generically. This second hypothesis 
could be seen as an extension of the Bohigas-Giannoni-Schmit conjecture [33] to the 
statistical properties of the spectrum of the classical evolution operator. This way, 
the random matrix properties in fully chaotic systems would have a fully classical 
counterpart. 

In semiclassical theories, the quantum correlations can be related to correlations 
acting among the actions of the periodic orbits [2] . In scaling systems like billiards, 
the action coincides, up to a constant factor, with the period (or length). There- 
fore action correlations are equivalent to period correlations. Since, as we have here 
shown, the density of periods of periodic orbits can be expressed in terms of the 
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eigenvalues of the classical evolution operator, it follows that the correlations be- 
tween periods of periodic orbits can be expressed in terms of correlations acting 
among the Ruelle-Pollicott resonances. The quantum spectral correlations are thus 
mapped, via semiclassics, into classical spectral correlations. Through this con- 
nection, the RMT conjecture of the quantum fluctuations should have a classical 
counterpart, that applies to the fluctuation properties of the position in the com- 
plex plane of the Ruelle-Pollicott resonances. This gives some support to our second 
hypothesis. 

The two previous conjectures determine the gross features of the low lying spec- 
trum of the evolution operator in fully chaotic billiards and, as a consequence, of the 
long time behavior of the density and correlations of periodic orbits. The random 
matrix universality observed in quantum systems may have, by semiclassical argu- 
ments, a classical counterpart. The "ergodic" zero located at the origin certainly 
plays an important role [12]. We are here suggesting a clear and explicit additional 
link between the statistical properties of the classical and quantum spectrums, now 
involving the complex resonances. 

To conclude, we briefly discuss the spectrum of the evolution operator for discrete 
maps, and show that important qualitative differences with respect to smooth flows 
occur. We have in mind area-preserving classically chaotic maps acting on two- 
dimensional phase spaces, like for example the kicked Harper or the kicked top. The 
time now takes only discrete values, r = n, n = 1, 2, 3, . . . (in some arbitrary units), 
and the set of possible lengths of periods of periodic orbits is trivial, just integers. 
The "return probability" or trace of the evolution operator is still expressed as a 
sum over the periodic points [4], 

oo 

Rin) = tr/:- = J2Y. |det(M._i)| ^-". = E ^^ ^"^" ' (50) 

p r=l ' ^ P ' ' ^ 

where Up is the (integer) period of the periodic orbit p, and the 7's are the Ruelle- 
Pollicott resonances. The latter are the zeros of (PT|) . making the replacement Tp — > 
Up. Inverting Eq. (j5n|l as in 21 3. formula follows for the number of primitive periodic 
orbits of period n, 

Afin) = ^J2fi{m) |det(M„/„-l)| J^^?, e"^"/"^ . (51) 

m/n 7 

As in ^ we have ignored the fluctuations of the factor | det(M^ ~ 1)1 which may 
occur at short periods. An alternative formula, equivalent to (|T6|l . is obtained 
from the corresponding Ruelle zeta, A/'(n) = n~^^^,^{m) ^^g-q e''"/"^, where 
the ?7's are the poles and zeros of Eq.Q with the replacement mentioned above. 
Since the time t = n changes by unit steps, the smallest scale over which temporal 
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variations can occur is one. It follows from Eq.(|5(J|) that the complex resonances 
satisfy —tt < Im(7) < vr. Therefore, the natural variable to analyze discrete maps is 
not s, but rather z = e~^. By this transformation, the ergodic zero 70 = is located 
at z = 1, and other real and complex Ruelle-Pollicott resonances with Re(7) > 
lie inside the unit disk \z\ < 1. 

Substantial differences are expected between the structure of the spectrum of 
Z{s) and Zpi{s) for chaotic maps with respect to that of chaotic continuous flows, 
and in particular with respect to chaotic billiards. The main reason for that is 
the simplicity of the spectrum of periods of periodic orbits in the case of maps. 
Since that spectrum is made of integers, the only non-trivial information carried by 
Eqs. (l5(J|) and (j51|) concerns variations in the stability factors and number of orbits 
of a given period (e.g., average coarse-grained properties). This is in contrast with 
continuous chaotic billiards, where a much more subtle and delicate information 
is encoded in the spectrum of resonances, namely the nontrivial distribution of 
the periods of the orbits. In maps, that distribution collapses to a simple and 
highly degenerate spectrum. The simplicity of the average properties encoded in 
Eq. (j5(J|) . without fine-grained structure, implies a simpler spectrum of resonances. 
In particular, no concentration of resonances over a "critical" line is expected to 
occur (this would be a critical circle inside the unit disk in the z variable), since its 
presence is associated to a harmonic decomposition of the distribution of periods 
in -R(t), whereas R{n) strictly tends to a constant for long times. In hyperbolic 
maps, isolated resonances, without any special structure in the radial direction, are 
therefore expected generically inside the unit circle. No particular connection with 
random matrix theory is therefore established concerning the statistical properties 
of the resonances of chaotic maps. This picture seems to be confirmed by recent 
numerical simulations [14]. 

Semiclassically, the reason for the important differences with respect to chaotic 
billiards is also clear. Random matrix requires correlations between actions of peri- 
odic orbits. In scaling systems like billiards, the action of an orbit is proportional to 
its length (or period). Therefore action correlations translate into length (or period) 
correlations, which in turn, through Eq. lj^Hj) . translate into resonance correlations. 
In maps, the set of possible periods (which is trivial) is not related to that of actions 
(which is nontrivial). Therefore, the spectrum of the evolution operator, which is 
closely related to the spectrum of periods of the periodic orbits, looses its connection 
with random matrix theory. 

References 

[1] W. Parry and M. Pollicott, Ann. of Math. 118 (1983) 573. 



26 



[2] 
[3] 
[4 

[5] 
[6] 
[7 



[9 

[lo; 
[11 

[12; 
[is; 

[14 

[is; 

[16 

[1?: 



N. Argaman, F.-M. Dittes, E. Doron, J. P. Keating, A. Yu. Kitaev, M. Sieber, 
and U. Smilansky, Phys. Rev. Lett. 71 (1993) 4326. 

M. PoUicott, Invent. Math. 81 (1985) 413; D. Ruelle, Phys. Rev. Lett. 56 (1986) 
405. 

P. Cvitanovic, R. Artuso, R. Mainieri, G. Tanner and G. Vattay, Classical and 
Quantum Chaos, www.nbi.dk/ChaosBook/, Niels Bohr Institute (Copenhagen 
2001) 

P. Gaspard, Chaos, Scattering and Statistical Mechanics, Cambridge University 
Press, Cambridge UK, 1998. 

M. C. Gutzwiller, Chaos in Classical and Quantum Mechanics, Springer Verlag, 
New York, 1990. 

R. Bahan and C. Bloch, Ann. Phys. (N.Y.) 69 (1972) 76. 

His remarks makes reference to the semiclassical connection between the set of 
periodic orbits and the quantum eigenvalues [23]. 

P. Dahlqvist, Physica D 83 (1995) 124. 

O. Agam, B. L. Altshuler and A. V. Andreev, Phys. Rev. Lett. 75 (1995) 4389. 

K. Pance, W. Lu and S. Sridhar, Phys. Rev. Lett. 85 (2000) 2737; S. Sridhar 
and W. T. Lu, J. Stat. Phys. 108 (2002) 755. 

A. V. Andreev, O. Agam, B. D. Simons and B. L. Altshuler, Phys. Rev. Lett. 
76 (1996) 3947. 

O. Bohigas, P. Leboeuf and M. J. Sanchez. Found, of Phys. 31 (2001) 489. 

J. Weber, F. Haake, P. A. Braun, C. Manderfeld and P. Seba, J. Phys. A 34 
(2001) 7195; S. Nonnenm acher, arXiv :nlin. CD/03 01014, I. Garcia-Mata, M. 



Saraceno and M. E. Spina, arXiv:nlin.CD/0301025l 



R. Aurich and F. Steiner, Phys. Rev. A 45 (1992) 583. 

S. Smale, The Mathematics of Time, Springer- Verlag, New York, 1986. 

D. Ruelle, Statistical Mechanics, Thermodynamic Formalism, Addison- Wesley, 
Reading, MA, 1978. 



[18] H. M. Edwards, Riemann's Zeta Function, Academic Press, New York, 1974. 



27 



[19] N.X. Chen, Phys. Rev. Lett. 64 (1990) 1193; J. Maddox, Nature (london) 344 
(1990) 377; D. Ming, T. Wen, J. Dai, X. Dai, and W. E. Evenson, Phys. Rev. 
E 62, (2000) R3019; N. V. Surovtsev, Phys. Rev. E 64, (2001) 061102. 



[2o; 

[21 

[22 

[23; 

[24 
[25; 

[26; 



[27 
[28; 
[29 

[3o; 

[31 

[32; 

[33 



P. Cvitanovic and B. Eckhardt, J. Phys. A 24 (1991) L237. 

H. H. Rugh, Nonhnearity 5 (1992) 1237. 

B. Eckhardt, in Quantum Chaos (Eds. G. Casati, I. Guarneri, U. Smilansky), 
CXIX International School "Enrico Fermi", North-Holland, Amsterdam, 1993, 
p.77. 

A. M. Ozorio de Almeida, Haniiltonian Systems: Chaos and Quantization, 
Cambridge University Press, Cambridge, 1988. 

J. Hannay and A. M. Ozorio de Almeida, J. Phys. A 17 (1984) 3429. 

Ya. V. Pesin, Dokl. Akad. Nauk SSR 226 (1976) 774 (Soviet Math. Dokl. 17 
(1976) 196. 

If the leading resonances would have satisfied g^ > htop, then the leading r -^ oo 
asymptotic oscillatory contributions to the density would have been terms with 
the largest possible value of m in the sum. But this is inconsistent with the fact 
that repetitions do not dominate the proliferation of periodic orbits. 

N. L. Balazs and A. Voros, Phys. Rep. 143 (1986) 109. 

M. C. Gutzwiller, Phys. Rev. Lett. 45 (1980) 150. 

A. Selberg, J. Indian Math. Soc. 20 (1956) 47-87; D. Hejhal, Duke Math. J. 43 
(1976) 441-482; J. Bolte, Open Sys. and Information Dyn. 6 (1999) 197. 

D. Ruelle, in Mathematical Physics X (ed. K. Schmiidgen), Springer- Verlag, 
Berlin, 1992, p.43. 

M. V. Berry, in Quantum Chaos and Statistical Nuclear Physics (Eds. T. H. 
Seligman and H. Nishioka), Lectures Notes in Physics 263, Springer Verlag, 
Berlin, 1986, p.l. 

A. Connes, Sel. Math., New Ser. 5 (1999) 29. 

O. Bohigas, M.-J. Giannoni, and C. Schmit, Phys. Rev. Lett. 52 (1984) 1. 



28 



